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ABSTRACT 


The  plane  metal  sla^  problem  hae  been  treated  in  a  rigoroue  manner 
for  constant  heat  input  and  the  4th  degree  nonlinear  radiation  boundary 
condition  at  the  front  surface  and  an  insulated  rear  surface.  This  has  been 
possible  by  formulating  the  problem  as  a  nonlinear  Volterra  integral 
equation.  An  approximate  solution  has  been  constructed  with  appropriate 
error  bounds.  This  technique  shows  that  the  solution  with  the  corresponding 
linear  boundary  condition  yields  an  accurate  solution  for  a  large  number  of 
problems  of  practical  interest.  Numerical  examples  including  rigorous 
estimates  of  their  validity  are  given  for  several  cases.  For  the  case  where 
this  approach  does  not  yield  sufficiently  accurate  results*  two  alternate 
approaches  are  suggested. 
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HEAT  TRANSFER  TO  A  METAL  SLAB  WITH 
RADIATION  BOUNDARY  CONDITIONS 


PURPOSE 

The  problem  of  constant  heat  input  to  a  planar  metal  slab  with  the 
4th  degree  nonlinear  radiation  boundary  condition  at  the  input  surface  and 
an  insulated  rear  wall  is  considered  herein.  The  purpose  of  this  paper  is 
to  develop  an  approximate  solution  from  the  rigorous  formulation  of  the 
problem  and  to  demonstrate  its  range  of  validity  with  appropriate  error 
bounds.  This  is  accomplished  by  formulating  the  problem  as  a  nonlinear 
Volterra  integral  equation  and  working  with  the  classical  iterated 
solutions.  The  properties  of  the  solution  are  developed  from  both  a 
physical  and  a  mathematical  approach. 


INTRODUCTION 

In  this  section  the  problem  will  be  stated  in  the  customary  form  of 
a  partial  differential  equation  with  boundary  conditions;  then  physical 
arguments  for  the  behavior  of  the  solution  will  be  presented.  The 
arguments  concerning  the  behavior  are  confirmed  in  Appendix  II,  where 
the  same  conclusions  are  drawn  from  the  integral  equation  formulation 
of  the  problem.  The  integral  equation  formulation  of  the  solution,  which 
has  the  advantages  of  including  the  boundary  conditions  directly,  is 
presented.  It  is  based  on  the  derivation  in  Appendix  I.  The  technique 
for  solving  the  integral  equation  is  shown  and  the  basic  approximations 
are  given. 

The  geometry  and  coordinate  system  of  the  problem  are  depicted 
in  Fig.  1.  The  differential  equation  and  boundary  conditions  are  given 
in  Eq.  (1). 


(1) 


- —  -  -  =  V  ;  T  ^  V  ,  V 

hi 

=  -  S2L  +  \  U*  (O.t) 

d| 
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whe  re 


aU(l>T)  =0 

hi 

u(e»0)  =  U(O.O) 

Qo  =  qo  +  £  «■  u^(o.o) 

Kt  ,  z 
T  =  - J-  ,  I  =  — 

I  i 

€  <ri 

is  the  heat  input 
e  is  the  emissivity 
(T  is  the  Stephans -Boltzman  constant 
t  is  time. 


Fig.  1.  Geometry  of  metal  problem,  where  k  is 

the  thermal  conductivity  and  K  the  diffusivity. 
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In  order  to  determine  the  behavior  of  the  temperature  when  subject  to  the 
nonlinear  boundary  condition,  the  boxmdary  condition  at  the  front  surface  is 
examined  in  detail.  In  Eq.  (2)  the  boxmdary  condition  is  presented  in  a 
rewritten  form: 

(2)  =  -  3aL  +\  [U*(0,t)-U*(0,0)]  , 

H  k 

where  the  term  X  [U*(0,  t  )-U*(0,0)]  represents  the  heat  loss  due  to  radiation. 

If  the  nonlinear  term  X  U*(0,  t  )  were  not  present,  one  would  expect,  since  the 
input  is  contant,  that  the  temperature  would  increase  monotonically  with 
increasing  t  and  become  unbounded  as  t—»oo.  However,  with  the  nonlinear 
term  present,  as  the  temperature  increases  the  heat  input  decreases  due  to 
the  radiation  losses,  and  the  input  term  thus  decreases.  Hence,  as  before, 
the  temperature  would  be  expected  to  rise  monotonically  but  with  decreasing 
slope  as  T  is  increased,  and  finally  as  t-*  »  the  input  term  on  the  right  of 
Eq.  (2)  would  approach  zero.  At  this  point  there  would  be  no  net  heat  transfer; 
that  is,  the  quantity  of  energy  supplied  to  the  slab  would  equal  the  quantity  of 
energy  lost  through  radiation.  Thus  an  equilibrium  temperature  would  be 
reached.  This  temperature  will  be  denoted  by  Ug  and  is  obtained  by  equating 
the  right  side  of  Eq.  (2)  to  zero: 

(3)  Ug  -  [Qg/ecr]*/'* 

Furthermore,  it  is  expected  that  for  |  ^  0  U(|,t  )  xvill  rise  monotonically  from 
U(0,0),  the  initial  temperature  (which  is  assumed  to  be  less  than  Ug),  and  will 
approach  Ug  as  t-^oo.  The  temperature  in  the  interior  (0  <  §<  1)  will  be  less 
than  the  surface  temperature. 


Since  none  of  the  standard  techniques  for  partial  differential  equations 
seem  applicable  to  the  solution  of  Eq.  (1),  an  attempt  is  made  to  formulate 
Eq.  (1)  as  an  integral  equation.  This  was  begun  by  postulating  a  solution  of 
the  form 

1  » 

G(|,t  Ii',  tM  U*(e',T')  , 

0 

where  G(|,  t  \  is  determined  by  the  condition  that  U{|,t)  must  satisfy 

Eq.  (1)  and  the  boundary  conditions.  The.  term,  Uo(|.t),  is  a  solution  to 
Eq.  (1)  with  the  nonlinear  boundary  condition  at  ^-0  replaced  by  a  linear 
boundary  condition,  which  can  be  handled  by  standard  techniques.  The 
boundary  condition  chosen  was  that  which  corresponded  to  the  same  problem 
but  with  radiation  neglected,  namely. 


(4) 


u{e.T)  =  u„(e,T)-x 


U. 

€  0  T* 


(5) 
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di  k 
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The  solution  for  Uq^, t  )  is  well  known,  ‘  and  is  given  by 


(6) 


Uo(I,t)  =  u(o.o)  +  ^  [t  +  {  j  -  e+L  I*  } 


2  2 


2  y  -2  fit  1 

— j  >  {n  cos  nir  §)  e  J 


n=l 


I  /  f 

The  function  G(^,  t  |  |  ,  t  )  was  determined  and  details  of  the  derivation 
are  given  in  Appendix  I. 


(7) 


[l +2V  (cosmrQe  ”  ^  ^]U^{0.tM 

t‘'=  0  1^1 


In  the  integrand  of  Eq.  (7)  only  the  surface  temperature  is  unknown, 
thus  if  the  surface  temperature  can  be  determined,  then  the  tempera¬ 
ture  fields  throughout  the  slab  can  be  determined.  An  integral  equation 
in  the  surface  temperature  is  now  obtained  from  Eq,  (7)  by  setting 
4=0.  The  resulting  expression  is  given  by 


(8) 


u(t  ) 


g(T 


-  t'  )  u*  (t') 


where 


u(t  )  =  U{0,T  ) 
u„(r)=U<,(0.T)^ 


g(T  -r')  =  1  +  2 


I 


n-1 


T  -t'  >  0 


=  0  ,  T  -t'  <  0  . 

Before  proceeding  further,  the  relation  between  Uo(t  )  and  g(T  -t')  will 
be  given  in  order  that  Eq.  (8)  may  be  recast  in  more  illuminating 
form.  First  Uq(t  )  is  written  in  the  form 
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(9)  Uo(t  )  =  Uo(0)  +  -^  i);jT) 

00 

1  2  •Z  *r 

where  )  =  t  +  —  -  — —  )  n  e’  , 

O  3  IT*  Zj 

n=l 

and  although  g(T  -t^)  is  not  bound,  it  is  integrable  and  an  integration 
shows 


(10)  g(T-T') 

=0 

Using  this  relation,  Eq.  (8)  is  now  recast  in  the  form 

T 

(11)  u(t)  =  Uo(0)  +  \  J  g(T-TM  [u*  -  U^t')] 

t'  =0 

where  u^  is  given  by  Eq.  (3).  The  details  of  this  derivation  are 
given  in  Appendix  1.  The  properties  of  u(t  )  are  derived  in  Appendix 
II,  and  the  conclusions  of  this  work  show  that  for  Uq(0)  <  u^  , 


(12)  u(t  )  is  monotone  increasing 

«o(0)  <  u(t  )  <  Ug 
u(t  ) - S-Ug 


These  properties  are  shown  pictorially  in  Fig.  2. 


Equation  (11)  is  solved  by  an  iterative  technique  given  by 


(13a)  ui  (t  )  =  u^(0)  +  \J  g(T  -t'  )  [u/  -  u^*(0)] 

t'=0 

T 

(13b)  Un(T  )  =  Uo(0)  +  \  J  g(T  -t")  [ug*  -  uJ_j  (t')] 

t'=0 
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Fig.  2.  Pictorial  representation  of  u(t  )>  the  surface 
temperature  at  ^  '  0,  the  input  surface. 

where  it  can  be  shown  Uj^(t  ) - ►u(t  )  under  suitable  restrictions.  These 

are  discussed  briefly  in  Appendix  III»  but  the  chief  concern  in  this  report 
is  ui  (t  )  since  it  can  be  obtained  explicitly  and  because  of  its  relationship 
to  Uq(t  ).  The  computation  of  Ujj(t  )  for  n  ^  2  cannot  be  carried  out 
explicitly,  and  numerical  methods  must  be  employed. 

However,  as  suggested,  uj  (t  )  is  expressable  rather  simply, 
since  u^^  -  (0)  is  a  constant.  This  can  be  seen  from  Eq.  (10). 

Hence,  it  is  found  that  Uj  (t  )  is  given  by  the  expression 


(14)  u,  (t  )  =  Uq(0)  +  \  [u^*  -  Uq*(0)]  i|<^(t  ) 

The  expression  for  ui  (t  )  may  also  be  written  in  the  form 
(since  \  u^^  =  QqI  /k) 

(15)  ui  (t  )  =  u^(t  )  -  \  u^*(0)4»^(t  ) 


where  Uq(t  )  is  given  by  Eq.  (9)  and  is  the  solution  to  Eq.  (1)  with  the 
nonlinear  boundary  condition  omitted  as  described  by  Eq.  (5).  From 
Eq.  (14)  it  is  seen  that  ui  (t  )  is  monotone  increasing  and  unboxmded  as 
T  -^oo;  thus  it  will  not  be  a  useful  approximation  for  all  r  .  However, 
for  T  small  (the  precise  meaning  will  be  apparent  later),  ui  (t  )  will  be 
a  useful  approximation.  In  Appendix  U  it  is  shown  that  Uj  (t  )  satisfies 
the  following  inequalities: 
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(16) 


0  <  Ui  (t  )  -  u{t  )<  \  [  Uq*(t  )  -  u^*(0)]  )  . 


In  terms  of  Uq(t  ),  the  solution  with  radiation  neglected^  this  inequality 
reads «  where  u^  (t  )  as  given  in  £q.  (15)  has  been  used. 

(17)  X.  Uq^(O)  4'q(t)  <  Uq(t)  -  u(t  )  <  \  [ujj*(t  )]  v1<q(t). 

The  inequality  in  Eq.  (17)  indicates  how  close  Uq(t  )  is  to  u(t  )  and  will 
be  illustrated  for  several  examples  in  the  next  section.  Thus  the 
inequalities  in  Eqs.  (16)  and  (17)  define  the  range  of  r  for  which  the 
approximation  is  accurate,  and  the  degree  of  accuracy.  Numerical 
techniques  for  evaluating  i)<o(t  )  and  thus  Uq(t  )  are  given  in  Appendix 
IV.  Also  as  a  result  of  Eq.  (17)  it  is  seen  that  by  a  slight  rearrange¬ 
ment  one  obtains  an  upper  and  lower  bound  on  u(t  ),  the  exact  solution. 
Explicitly  this  is  given  by  the  expression 

(18)  Uo(t  )  -  X.  Uq^  (t  )  i1Io(t  )  <  u(t  )  <  Uq(t  )  -  \  u^^  (0)  4Jq(t  )  . 


These  upper  and  lower  bounds  on  u(t  )  are  shown  graphically  for 
several  examples  in  the  next  section.  In  the  examples  considered 
there  the  quantity  \  Uq^(O)  ijt  (t)  on  the  right  of  Eq.  (18)  was  quite 
small  and  hence  was  neglected.  Note  that  u^(t  )  is  still  an  upper 
bound  on  u(t  )  when  \  Uq*(0)  v|<q(t  )  is  replaced  by  zero. 


NUMERICAL  RESULTS  AND  EXAMPLES 

In  this  section  numerical  results  are  given  for  a  1/2 -inch  copper 
slab  and  a  1/2 -inch  stainless  steel  slab.  The  results  concerning  the 
accuracy  of  u^(t  )  as  an  approximation  to  u(t  )  are  summarized  in 
Table  I  when  u(t  )  is  near  the  melting  temperature  of  the  slab.  The 
errors  given  in  the  table  are  the  largest,  since  the  approximation 
improves  as  t  decreases.  Also  it  is  felt  that  there  is  no  need  to  study 
the  solution  u(t  )  when  u(t  )  exceeds  the  melting  temperature  of  the 
slab,  for  surely  the  assumptions  that  the  thermal  parameters  are 
independent  of  temperature  are  no  longer  valid  in  this  temperature 
range.  The  thermal  parameters  and  other  data  relating  to  the  two 
slabs  are  given  in  Table  II. 
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TABLE  I 


1/2- 

inch  Copper  Slab 

1/2 -inch  Stainless  Steel  Slab 

Melting  Temperature  1325  K 

Melting  Temperature  1800°K 

Max. 

Max.  Percent 

Max. 

Max.  Percent 

€ 

Error 

Error 

Error 

Error 

0.  03 

0.  02 

<  0.  ooi^i 

0.  03 

10*K 

0.  6% 

1.0 

175^K 

<3.8% 

0.1 

Sl^'K 

2% 

1.  0 

310°K 

21% 

Table  I  describes  the  accuracy  of  the  approximate  solution  near 
the  melting  point  of  the  slab. 


TABLE  II 

I  =  1.  27  •  10”*meters  Uo(0)  =  300°K 


Parameter 

Copper  Slab 

Stainless  Steel  Slab 

watts 

^  meter('^Kj 
(meter)* 

{''K)  (sec) 

watts 

(meter)* 

389 

1. 14  •  10"* 

2.  8  •  10*‘ 

14 

3.  87  •  10 

2.8  •  10‘ 

-*6 

0.  03 

1.  0 

0.  03 

0.  1 

1.  6 

X  (°K)"^ 

5.  62  '  10’** 

1.  54.10"“ 

1.  56*  10"‘* 

5.  21  •  10"** 

5.  21-  10"“ 

u,('»K) 

6500 

2840 

6350 

4700 

2640 

Table II  gives  the  parameters  associated  with  the  two  klabs. 
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In  Figs.  3  and  4  the  results  for  the  copper  slab  are  given,  but  only  for 
the  case  c  =  1.  0.  The  results  for  e  =  0.  03  are  very  good,  as  can  be 
seen  in  Table  I ,  and  no  attempt  is  made  to  graph  them.  The  worst 
error  is  only  0.  02°K  and  occurs  when  Uq(t  1325°K.  The  results  for 
the  stainless  steel  slab  are  given  in  Figs.  5,6,  and  7.  The  worst  error 
occurs  again  for  e  =  1.  0. 


1500 
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y  900 
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3001 
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Point  Of  C 
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/ / 

- 

U 

// 

// 

// 

f 

7 

/ 

/ 

_ 

8  12  16 
t( Seconds) 


Fig.  3.  Figure  3  shows  Uq(t)  for  the  1/2-inch  copper  slab.  For  €=1.0 

the  exact  solution,  u(t)  lies  above  the  dashed  line  and  below  Uo(t). 
For  €  =  0.  03  the  dashed  and  solid  curves  are  too  close  to  depict 
separately. 
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0  4  8  12  16 


t( Seconds) 

Fig.  4.  Figure  4  shows  the  error  for  the  1 /2-inch  copper  slab  when 

e  -  1.0.  This  is  the  worst  case;  that  i8>  for  all  other  parameters 
fixed)  as  e  increases  the  approximation  becomes  less  accurate. 


Three  suggestions  for  improving  these  results  will  now  be  given. 
Two  of  the  techniques  will  depend  on  a  numerical  integration  and  the 
third  on  solving  a  linear  problem  which  is  closer  to  the  actual  physical 
problem. 


In  deducing  the  upper  bovmds  given  in  Eqs.  (16),  (17)  and  (18) 
relatively  crude  techniques  were  used,  for  example,  since  i|<q(t  )  is 
monotone  increasing,  the  rather  crude  result  follows 
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2200 


t  (Seconds) 


Fig.  5.  Figure  5  shows  U(,(t)  for  the  1/2-inch  stainless  steel  slab.  For 
€  s  0.  1  the  exact  solution  lies  between  Uq(t)  and  the  dashed  curve 
labeled  0. 1.  For  e  =  1.  0  the  exact  solution  lies  between  Uo(t) 
and  the  dashed  curve  labeled  (  -  1.0.  For  c  =  0.  03  the  dashed 
curve  is  too  close  to  Uq(t)  to  depict  seperately  on  this  graph. 

Hence,  it  is  quite  possible  that  the  large  discrepancy  indicated  for  the  stain¬ 
less  steel  example  is  an  unnecessarily  harsh  evaluation.  Therefore,  the 
first  approach  would  be  to  integrate  numerically  the  right  hand  side  of  the 
inequality  given  in  £q.  (19)  below. 


(19)  0<  u,(t)  -  u(t)<  X  j  g(T-T')  [uNt')  -  u^*(0)] 

t'=0 

This  result  is  deduced  from  £q.  (2-6b)  and  the  fact  that  u(t  )  <  Uq(t  )  . 
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0  4  8  12  16 

t( Seconds) 

Fig.  6.  Figure  6  shows  the  error  for  the  1 /2-inch  stainless  steel  slab  for 

€  =  0.  03  and  €  =  0.  1.  The  error  for  6  s  1. 0  is  obtained  by  multiplying 
the  ordinate  by  the  factor  10.  Thus  again  it  is  seen  that  as  €  increases^ 
so  does  the  error. 

If  this  approach  fails  to  give  the  desired  accuracy,  the  next  step  would 
be  to  compute  several  of  the  iterated  solutions  Uj^(t  )  for  n  >  2  and,  in  this 
exaiiiple,  for  0  <  t  <  4  seconds.  This  solution  for  0  <  t  <  4  seconds  is  used 
to  continue  the  solution  for  t  >  4  seconds.  This  is  accomplished  by  rewriting 
£q.  (8)  as  follows 

T  (t=4  sec)  T 

u(t  )  =  Uo(t  ) -xj'  g(T-T')  U*(T')-x/'g(T-T')  U*(t') 
t'=0  t'=t  (t=4  sec) 

The  first  two  terms  on  the  right  hand  side  of  the  equality  are  known. 

Considering  these  two  terms  as  the  initial  approximation  the  temperature 
for  t  >  4  seconds  can  now  be  computed  by  iteration.  However,  since  this 
approach  depends  quite  heavily  on  numerical  integration  the  resultant  errors 
from  the  integration  may  make  this  approach  prohibitive. 
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Fig.  7.  Figure  7  shows  the  maximum  percent  error 
for  the  1/2-inch  stainless  steel  slab  when 
€  =  1.0,  the  worst  case. 


The  third  approach,  as  mentioned  earlier,  consists  of  solving  a  linear 
problem  which  is  closer  to  the  physical  problem  and  evaluating  its  accuracy 
by  the  techniques  presented  here.  The  first  step  is  to  expand  the  right  hand 
side  of  Eq.  (2)  in  a  power  series  about  U(0,0).  This  is  given  by 

(20)  ^  (0,t)  =  -  -52i  +  4  XU*(0,0)[U-U(0,0)]  +\  H(U) 

ae  k 

where  H(U)  =  6  U*  (O,  0)[U-U(0, 0) ]*  +  4U(0, 0)[U-U(0, 0)] »+ [U-U(0, 0)]* . 

It  is  seen  that  the  right  hand  side  of  Eq.  (20)  involves  terms  which  are  linear 
in  U  and  terms  which  are  nonlinear  in  U.  The  nonlinear  terms  are  grouped 
in  the  factor  H(U).  It  is  possible  to  obtain  a  useful  solution  to  Eq.  (1)  when 
the  nonlinear  terms  are  neglected.  This  solution  has  two  advan';ages  over 
linear  solution  which  has  been  used  in  the  examples  considered  here. 
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(1)  it  take  I  into  account  some  radiation  loss  for  t  small 
while  the  previous  did  not. 

(2)  this  solution  will  be  bounded  as  while  the 

previous  was  not. 

Also,  the  nonlinear  term  may  be  expanded  in  a  power  series  about  the 
equilibrium  temperature.  This  will  yield  an  approximate  solution  which 
will  be  useful  for  t  large. 

No  numerical  results  based  on  this  approach  are  available  but  of  the 
three  methods  suggested  here  this  appears  to  be  the  most  promising. 


CONCLUSIONS 

The  metal  slab  problem  with  the  nonlinear  boundary  conditions  has 
been  formulated  as  a  nonlinear  Volterra  integral  equation.  A  simple  approxi¬ 
mate  solution  has  been  found  which  is  shoum  to  be  quite  adequate  for  several 
cases  of  practical  interest.  Error  bounds  indicating  the  validity  of  the 
approximation  have  been  given.  Techniques  have  been  suggested  for 
handling  the  case  when  the  error  becomes  too  large. 

Since  it  has  been  possible  to  determine  the  equilibrium  temperature 
by  examining  the  integral  equation  satisfied  by  the  surface  temperature, 
the  possibility  exists  that  the  equilibrium  temperature  for  more  complex 
shapes  can  be  determined  by  a  similar  procedure. 
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APPENDIX  I 


In  this  section  an  argument  is  given  for  reformulation  Eq.  (1)  as  an 
integral  equation.  A  similar  argument  based  on  Green's  second  theorem 
can  be  given.  The  equations  to  which  a  solution  is  sought  are  restated  in 
Eq.  (Al-la)  to  Eq.  (Al-ld). 

(Al-la)  .LJL  -  ^=0,  T>0,0<e<l 

S  I 


(Al-lb)  iH.(0,T)  =  -  Htti  +x  U^(0,t), 

ae  k 


(AI-lc)  ^(1,T)  =  0, 


I 

I 

I 


(Al-ld)  U(g,0)  =  U(0,0)  , 

where  the  symbols  are  defined  in  Eq.  (1).  An  attempt  is  made  to  seek  an 
integral  equation  formulation  of  the  problem  by  asking  for  a  representation 
of  U(|,  T  )  in  the  form 

so  1 

(Ai-2)  u(e,T)=Uo(e,T) -xj  J  G(e,T  ir.T')  u«(e'.T') . 

■K-o  r=o 


where  Uo(6»  t  )  satisfies  Eq.  (AI-1)  with  \  =••  0,  andG(|.T  t^)  is 
determined  by  the  requirement  that  U(^t)  satisfy  Eqs.  (AI-1).  The 
solution  for  t  )  is  well  known*  and  is  given  in  Eq.  (6). 

Th  requirement  that  U(^t)  satisfy  Eq.  (Al-la)  requires  that 
G(6,t  k'.r')  satisfy 


T"  (fc  T  ~  ($.  T  I t'  )  =  0 

;  OT 


T  >  0  0  <  Cl  1 

t'>0  oiC'l  1 
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Notice  that  satisfies  Eq.  (Al-la).  Saving  the  boundary  condition  in 

Eq.  (Al-lb)  until  last,  the  boundary  condition  in  Eq.  (AI-lc)  is  considered 
next.  First  note  that 


aup 


(l.T)  =  0 


» 


and  hence,  in  order  that  U(4,t)  satisfy  Eq.  (AI-lc),  it  is  required  that 


(i,T  le^.T^)  =  0  . 

Turning  next  to  the  initial  condition  in  Eq.  (Al-ld),  the  following 
conclusions  are  drawn.  Since 

Uo(e,0)  =  U(O.O)  , 

it  is  required  that 

G(S,0U'/)=0  . 

Now  attention  is  turned  to  the  nonlinear  boundary  condition  esqpressed  in 
Eq.  (Al-lb).  Next,  it  is  noted  that 


Qof 


—12.  (0,t)  =  - 

H  k 


Hence,  it  is  seen  that  if  G(&  t  \^,r  )  satisfies 


-  (0.T  tO  =  -5(6')  6(t-t')  , 


where  5(x)  is  the  Dirac  delta  function,  then  U(6»  t  )  satisfies  the  boundary 
condition  at  the  input  surface.  Here  use  has  been  made  of  the  fact  that 


00  1 


U^0,t)=J'  J  6(6')  6(t-t')  UVe',  t')  . 

/=0  6'=0 
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In  summary,  then,  it  is  seen  that  G(^t  satisfies  the  following 

system  of  equations: 


(AI-3)  -  ^=  0  0<g<l,  0<T 

S  ^  S  T  ” 

TT-  (0.  'k',TO  =  -5(|')6(T-T') 


—  (l.T  li'  t')  =0  ,  T>0 

G(e.ok'.  T^)  =  0 ,  o<i<i , 

for  ^5,5^5.^  0<T^  . 


The  solution  to  this  system  (Eq.  (AI-3))  is  quite  similar  to  that  for  Uo(^t), 
as  can  be  seen  by  comparison.  The  solution  for  G(^  t  k^,  t^)  is  given  by 

1  +  2  >  (cos  mr^Je  '  ' 

n=l 


(AI-4)  G(S,Tk'.T')  =  6(e') 


s{t  -t^) 


where 


e{r-r')  =  0,  t-t'  <  0 
=  1  ,  T  -T^  >  0 


The  formal  solution  for  G(^  t  obtained  by  the  Laplace  Transform 

techniques.  Substituting  this  expression  (Eq.  (AI-4))  into  Eq.  (AI-2)  and 
carrying  out  the  integration  with  respect  to  g  ,  one  obtains  the  result  in 
Eq.  (7)  of  Section  B,  namely, 


(AI-5) 


U(S.t) 


=Uo(i,T)-xJ 


t'=0 


SO 

l+zV  (cos  nv  Q  ^ 


n=l 


g(& 


T  -t'  )  =  1  +  2  ^  cos  niT  I  e"”  > 

n=l 


--  0 


T  -T^  <0  . 
17 
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It  has  not  been  possible  to  show  in  a  rigorous  fashion  that  U(^t  )  is 
determined  by  Eq.  (AI-5);  however,  the  argument  appears  reasonable, 
and  as  will  be  seen  in  the  next  section,  the  temperature,  U(4»t  )^has 
all  of  the  expected  physical  properties. 
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APPENDIX  II 


In  this  section  the  properties  of  u(t),  the  surface  temperature,  are 
established.  These  properties  are  deduced  by  assuming  that  u(t)  is 
bounded  and  integrable  and  Uo(0)<  u®.  The  case  Uo(0)  =  u^  has  the  trivial 
solution  u(t)  a  Ug  and  the  case  u^CO)  >  Ug  is  of  no  interest  physically, 
although  it  can  be  handled  in  the  same  manner  as  the  present  case. 

From  Eq,  (AI-5)  in  Appendix  I,  by  setting  5=0  one  obtains  the 
following  expression: 


(AlI-1) 


where 


“(t  )  =  Uq(t  )  -  g(T-T')  u^(t') 

r*  =0 


u(t)  =U(0.t) 


g(T-T0 


=  Ug(0,T) 


=  1+2 


^  2  2 

•n  it*(t  -t^) 
c 


r 

n=l 


=  0 


T -t'>  0 
T -T^  <  0  , 


Note  that  Eq.  (AII-1)  implies  that  Ug(T )  >  u(t  ),  t  >  0,  ElemenUry 
theorems  show  that  g(T  -t')  is  summable  and  that 

T 

(AII-2)  4'o(t)  =  y'  g(T-T')  . 

t'=0 


With  this  information  and  the  relation  between  (t  )  and  u  (t  ),  one 
deduces  the  integral  equation  given  in  Eq,  (11),  ° 


(AII-3)  u(t)  =  Uo(0)  +\J*  g(T-T')[Ug*  -  u*(t')] 

t'  =0 


where 


u 


4  - 


Qo  QqI 
€  O’  "  X  k 
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Several  elementary  observations  are  now  made. 


u(0)  =  Uo(0), 

and  there  exists  a  t*  such  that  u{t  )  is  continuous  for  0<  t  <  t*  .  Since 
u(t  )  is  continuous  for  t  <  t  *  and  u(0)  =  Uo(0)  <  u^,  it  follows  that  u(t  )  <  Ug 
for  some  interval  to  the  right  of  t  ~  0.  Thus,  it  also  follows,  since  the 
integrand  in  Eq.  (AIl-3)  is  positive,  that  u(t  )  is  monotone  increasing  in  the 
same  interval  that  for  which  u(t  )  <  Ug.  Thus  it  has  been  established  that 


u(t  )  is  monotone  increasing  t  <  t  o 
where  u(t  q)  - 

It  will  now  be  established  that  t  ^  is  not  finite.  For  t  >  t  g, 
examination  of  Eq.  (All* 3)  shows  that 

To 


(An-4) 


u(t  )  =  Uo(0)  +  \J  g(T  -t'  )[Ug^ 


./ 


;  0 


u^tM] 


g(T-T')  [m* 


Since  t  >  t  q,  it  follows  that  g(TQ-T^)  >  g(T-T^),  Hence,  since  u(t')  <  Ug 

for  t'  <  T  o  , 


(AII-5) 


g(T-T' )lue'‘  "  "*(t')]»  t>t, 


But  it  is  assumed  that  u*(t')  >  Ug^  for  t'>  t  q,  which  implies  u(t  )  <  Ug  ; 
this  contradicts  the  assumption  that  u(t  )  equals  Ug  for  finite  t  .  Thus  u(t  ) 
is  a  monotone  increasing  function  for  t  ^  0,  bounded  below  by  Ug(0)  and 
above  by  Ug,  and  further 


u(t> 


u. 


These  are  precisely  the  results  based  on  physical  reasoning. 
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Next  the  bounds  on  the  difference,  (t  )  •  u(t),  will  be  derived. 
From  Eq.  (All*  3)  and  the  expression  for  (r )  given  in  Eq.  (AIl>6a),  the 
following  results  are  deduced! 


T 

{AII-6a)  Ui(t)  =  Uq(0)  +  X-  T  g(T-T')  [u^*  -  u^*(0)] 

t'=0 


Hence,  it  is  seen  that  (since  u(t  )  >  u^(0)) 


(AII-6b)  0<  u^(t)  -  u(t)  =  \j  g(T-TO  [u^(t')  -  Uq^  (0)1  . 

t'=0 


From  the  fact  that  u*(t  )  is  monotone  increasing,  and  bounded  by  the 
monotone  increasing  function  Uq^(t},  it  follows  that 


(An-7)  0  <  Uj  (t  )  -  u(t  )  <  X  [u^*(t  )  -  u^*(0)]  i|<^(t  )  . 


This  is  the  result  stated  in  Eq.  (16).  Thus  the  properties  of  u(t)  as  well 
as  the  relationship  between  u(t  ),  u^  (t  )  and  \i^{r  )  have  been  deduced. 
Some  of  the  properties  of  Un(T )  for  n  >  2  are  deduced  in  Appendix  UI, 
although  they  are  not  of  particular  interest  as  far  as  this  report  is 
concerned. 
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APPENDIX  m 


As  pointed  out  previously,  since  Uj  (t  )  is  unbounded  as  t— >oo,  it 
surely  will  not  offer  a  good  approximation  to  u(t  )  for  all  t  .  This  is  also 
illustrated  in  an  example  in  Section  2.  It  is  therefore  the  purpose  of  this 
section  to  show  how  the  iterated  solutions,  Un(T ) ,  n  >  2,  can  be  used  to 
improve  the  approximation.  In  this  section  it  will  be  shown  that  for  t 
less  than  a  certain  t  q,  the  convergence  of  Uj,(t  )  to  u(t  )  is  guaranteed. 
For  T  larger  than  this  value  another  method  is  suggested.  Let  Tj  be 
defined  as  follows: 

(AIII-l)  u,(tJ=u^  . 

The  next  iterated  solution  U2(t)  is  given  by  (as  defined  by  Eq.  (13b)) 

T 

(t)  =  Uq(0)  +  Xj'  g(T-T')  [U^^  -  Uj^(t')] 

0 


Then  for  t  <  t  i , 

Uo(0)  <  Uj(t  )  and  is  monotone  increasing. 
Further,  for  t  <  ti 

Uj  ( T  )  5.  '^1 

'^o(0)  1  ^  *  (t  )  <  u(t  )  <  u  j  (t  ) 


with  the  equality  signs  holding  for  x  =  0  only.  To  show  that  u(t  )  >  u,(t  ), 
consider  the  difference 


^(t)  -  Ua(T)  -  xT  g(T -t')  [ui*(t')  -  u*(t' )]  . 
t'=0 
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Since  Uj  (t)>  u(t),  t  <  the  result  follows.  In  a  similar  manner  one 
can  show 


u(t)  <  u^(t)  <  Uj(t),  t  <  Tj 


In  general  one  has  the  following  scheme  with  each  Uj^(t  )  a  monotone 
increasing  function  of  r  for  t<  T|  , 

I - ^ ^ ^ - hH - ^ ^ - 1 

'lotO)  u(t)  U7(t)  Ugir)  UjIt  )  \i^{r) 

Fig.  ni-1. 

In  this  figure  a  pictorial  representation  of  the  type  of  convergence  is  shown. 
The  question  to  ask  is,  for  what  range  of  t  does  the  sequence  |u(t  )-Un(T  )  | 
converge?  This  sequence  can  be  shown  to  converge  for  all  r  which  satisfy 


(Ain-2)  4X  4/q(t)  u,  *(t)  <  1 


Denote  the  largest  least  upper  bound  of  the  t 's  which  satisfy  Eq,  (AIII**2) 
by  Tj  :  then  for 

(Am- 3)  T  <  min  (tj  ,  T2  ]  the  sequence  ^u(t)  . 

n>4oo 

For  T  >  Tq  where  Tq  <  min  [tj  ,  tj  ],  another  technique  must  be  sought. 
Let  T  Q  be  a  value  for  which  an  accurate  iterated  solution  has  been  found. 
Rewrite  the  integral  equation  for  u(r )  as  follows: 

To  T 

u(t)  =  Uo(0)  +\j  g(T-T')  [u<^-U*(t')]  +\J  g(T-T')[Ug*-U^(T^)] 

=0  T '“T  _ 

O 

for  T  >  Tq. 

Note  that  the  term 
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t'  =0 


is  known,  and  hence  for  t  >  one  is  faced  with  solving  the  integral 
equation 


u(t)  =  f(T)  +xj  g{T-T')[Ug^“'l*(''’')l 


t'  =t. 


It  is  now  hoped  that  this  equation  can  be  iterated  to  obtain  approximate 
solutions  for  t  >  r^.  However,  this  matter  will  not  be  pursued  further 
here. 
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APPENDIX  IV 


In  this  section  some  information  concerning  the  calculation  of  i{>  (t  ) 
is  given.  This  function  is  given  by  ^ 

(AIV-l)  )  =  T  +  ^ 


and  due  to  the  exponential  terms  in  the  series  one  esqiects  good  convergence 
for  T  large.  This  is  indeed  the  case  as  the  following  tables  indicate. 
However,  for  t  small  the  convergence  is  slow  and  another  form  for  the 
expression  of  i|/  (t  )  is  sought.  This  esqpression  and  a  resulting  simple 
approximating  i&nction  is  found  by  noting  that  g{r )  is  aTheta  function. 
Relationships  between  theTheta  functions  yield  the  desired  relationship.’ 


(AIV-2)  g(T-T')  =0,(0.  liir(T-T')) 


where 


=  0,(o| - i - ) 

^  ir(T-T^  )■  ir  (t  -t'  ) 

0s  (0,  |ix)  = 

00 

1 + 2  y  '  * 

n=l 


The  following  results  were  then  deduced,  and  are  given  in  Tables  III,  IV, 
and  V.  Table  III  gives  the  general  results,  Table  IV  gives  results  when 
accuracy  to  five  places  is  needed,  and  Table  V  gives  results  only  when 
three-place  accuracy  is  needed. 
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TABLE  ni 


Various  approximate  expressions  for  )  and  bounds 
denoting  how  close  the  bounds  are. 


I 


1 

I 

1 
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I 

r 

L 

y 

L 

I 

1 

I 


TABLE  IV 


Gives  approximate  e3q>re88ion8  for 
accurate  to  five  significant  figures. 


Range  of  t 

Approximate  Expression  ) 

Five  significant  figures 

0  <  T  <  .1 

*<.<-)•  2  iT 

,  1  <  T  <  ,  15 

.  ,  ,  1  ^  ,  -tr*T  1  1  -9tt*t 

4/-(t  )  =  t  +  —  -  __  (e  +  —  e  +  T  *  ) 

^0  '  3  IT*  4  9 

.15  <  T  <  .  30 

.  /  »  -■n’*T  1  -4ir*T. 

*„(’■)“  T  +--^(«  +~  e  ) 

.30  <T  <1.25 

1  2  -Tr*T 

4'o(t)  =  t  —  (e  ) 

V  TT 

1.25<  T 

*^(t)=T  +1 

TABLE  V 

Gives  approximate  expressions  for  ) 
accurate  to  three  significant  figures. 


Range  of  t 

Approximate  Expression  ) 

Three  significant  figures 

0  <  T  <  .1 

.1  <  T  <  .2 

‘I'ot-r)  =  +  j)  -  -7*'  1 

.2  <  T  <  .7 

1  2  -Tr*T 

‘I'o(t  )  =  (t  -t)  -  — Y  ) 

TT 

.7  <  T 

)  =  T  +  i 

I 
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